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ABSTRACT 

We study the distribution of the star formation rate (SFR) and metallicity of damped 
Lyman-a absorbers (DLAs) in the redshift range z — — 4.5 using cosmological 
smoothed particle hydrodynamics (SPH) simulations of the A cold dark matter model. 
Our simulations include standard radiative cooling and heating with a uniform UV 
CS| ■ background, star formation, supernova (SN) feedback, as well as a phenomenological 

^ ' model for feedback by galactic winds. The latter allows us to examine, in particular, 

0^ , the effect of galactic outflows on the distribution of the SFR and metallicity of DLAs. 

We employ a "conservative entropy" formulation of SPH which alleviates numerical 
overcooling effects that affected earlier simulations. In addition, we utilise a series 
of simulations of varying boxsize and particle number to investigate the impact of 
numerical resolution on our results. 



o 

f-^ , We find that there is a positive correlation between the projected stellar mass 

■ density and the neutral hydrogen column density (Nhi) of DLAs for high Nhi systems, 

and that there is a good correspondence in the spatial distribution of stars and DLAs 
in the simulations. The evolution of typical star-to-gas mass ratios in DLAs can be 
characterised by an increase from about 2atz = 4.5to3atz = 3, tolOatz = l, and 
jr-H | finally to 20 at z = 0. We also find that the projected SFR density in DLAs follows 

the Kennicutt law well at all redshift s, and the simula ted values are consistent with 
the recent observational estimates of this quantity bv IWolfe. Prochaska. fc Gawiserl 
(|2003allbh . The rate of evolution in the mean metallicity of simulated DLAs as a 
function of redshift is mild, and is consistent with the rate estimated from observations. 
/\i 1 The predicted metallicity of DLAs is generally sub-solar in our simulations, and there is 

H , a significant scatter in the distribution of DLA metallicity for a given Nhi- However, 

we find that the median metallicity of simulated DLAs is close to that of Lyman- 
break galaxies, and is higher than the values typically observed for DLAs by nearly an 
order of magnitude. This discrepancy with observations could be due to an inadequate 
treatment of the SN feedback or the multiphase structure of the gas in our current 
simulations. Alternatively, the current observations might be missing the majority of 
the high metallicity DLAs due to selection effects. 

Key words: cosmology: theory - galaxies: evolution - galaxies: formation - methods: 
numerical. 



1 INTRODUCTION 

Damped Lyman-a absorbers (defined as quasar absorption 
systems with column density Nhi > 2 x 10 20 cm -2 ) have 
similar neutral hydrogen column density as galactic disks in 
the Local Universe. They are high concentrations of neutral 
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gas which could form st ars and eventually evolve into galax- 
ies that we see today (IWolfe et al.1 1986 V While the most 
numerous absorbers are those with low column densities, 
the high column densities of DLAs more than compensate 
for their relative paucity when the relative contribution to 
the total neutral hydrogen (Hi ) mass density is considered. 
It is thought that DLA s dominate the Hi gas content of 
the Universe at z ~ 3 (lLanzettal Il993l : IWolfe et alJll995l : 
IS torrie- Lombard! & Wolfe 2000), and thus they are serving 
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as an important reservoir of neutral gas for star formation at 
high redshift. Therefore, studying the physical properties of 
DLAs, such as the distribution of their star formation rates 
and metallicity, will provide us with important complemen- 
tary information to those provided by the emitted light from 
stars in high-redshift galaxies. 

It has become cle ar from recent obser vational (e.g. 
lAdelbererer et all Il998l: ISteidel et"aD Il99<t) and theo- 
retical (e.g. iMo fc Fukupdta "Tl99d: iBaugfi et al.1 Il998t 
Mo. Mao fc White! Il999t Katz. Hernauist fc Weinberg! 
19991: iKauffmann et all Il999l: iNagaim^e! |200j: 
Weinberg. Hernauist. &: Katzl [2002) studies of Lyman- 
break galaxies (LBGs) at z ~ 3 — 4, that the assembly 
of galaxies is actively going on at z ~ 3, consistent with 
what is expected from the hierarchical structure formation 
paradigm based on a cold dark matter (CDM) model. 
The star formation rates of LBGs at z = 3 — 4 are now 
commonly observed using emission lines such as Ha, or 
by fitting t he entire spectrum w ith population synthesis 
models (e.g. IShaolev et alll200lh . These measurements of 
SFRs typically give values of ~ 50 M©yr _1 for LBGs, 
allowing one to place a lower limit on the cosmic star 
formation rate by combining the information on the SFR 
and number density of LBGs. 

Furthermore, many theoretical and observational 
studies suggest that the cosmic star formation rate 
rises towards high redshift, even beyond z — 3 (e.g. 
Pa^c^eJte^Larrzje^^ 

19991: iNagamine et all l200ll: lLanzetta et all |2002t 
Soringel fc Hernauist! l2003bt iHernauist fc Springel! l200fl . 
Therefore, the conversion of gas into stars is taking place 
at a significant rate in these high-redshift galaxies at z > 3. 
Studying this conversion process can help us understand the 
galaxy formation better, and possibly constrain structure 
formation scenarios such as the hierarchical CDM models. 

However, estimating the star formation rate in DLAs 
has proven to be very difficult, because the stellar coun- 
terparts of D LAs cannot be detected in emission due to 
their faint ness. Wolfe, Prochaska, & Gawiser (2003a b) have 
therefore adopted a different approach and recently esti- 
mated the SFR in DLAs in the redshift range z — 3 — 4 
by utilising the Cn* 1335.7 absorption line. They infer the 
heating rate by equating it to the cooling rate measured 
from Cn* absorption. Since the heating rate is proportional 
to the SFR, they can then deduce the SFR from the esti- 
mated heating rate. The values of SFR they meas ure by this 
meth od roughly agree with the Kennicutt law (|Kennicuttl 
1998) with some scatter. The robustness of this technique 
still needs to be tested with future observations, but it is 
encouraging that observers are making progress towards a 
'direct' measurement of the SFR in DLAs. 

On the other hand, there is a wealth of observational 
data on the chemical abundances of DLAs from detailed 
studies of quasar absorption line spectra (see Pettini 2003 
for an excellent review). DLAs are generally inferred to be 
metal poor at all redshifts fe.g.lPettini et al| ll994l: iLu et all 
ll99d : IPettini et al.lll999l : IProchaska fc Wolfell99<ft . The col- 
umn density- weighted mean abundance of zinc is measured 
to be ~ 1/13 of solar (([Zn/H]) = —1.13), which suggests 
that DLA galaxies (i.e. galaxies in which DLAs arise) are 
at early stages of their chemical evolution. The measured 
values of [Zn/H] span more than an order of magnitude at 



each redshift, and the chemical enrichment appears to pro- 
ceed at different rates in different DLAs. The metallicity 
distribution of DLAs peaks at a value in-between those of 
halo and thick-disk stars of the Milky Way (Fig. 8 of Pettini 
2003; Wyse & Gilmore 1995; Laird et al. 1988), suggest- 
ing that the kinematics of DLAs should be closer to stellar 
haloes and bulges than to rotationally s upported disks. This 
findin g is at odds with the proposal bv lProchaska &; W olfe 
( 1998) that most DLAs are large disks with rotational veloc- 
ities in excess of 200 kms -1 . An interesting fact is that there 
is little evide nce for any redshift evolution in the metallic- 
ity o f DLAs (IPettini et al.lll999l : IProchaska fc WolfdExxl 
2002), despite the fact that most analytic models of global 
chemical evolution of the Universe predict substa ntial evo- 
lutio n in metallicity as a fun ction of time (e.g . |P ei &; Fal3 
Il995ll However, most recently IProchaska et al J ([2003F) have 
detected mild evolution with a significantly larger sample 
of DLAs. It is interesting to see if full cosmological simula- 
tions agree with observations or analytic models of global 
chemical evolution, and what they predict in detail for DLA 
metallicities. 

Theoretical understanding of the physical prop- 
erties of DLAs still remains primitive. Several 
semi- analytic models o f DLAs a nd th ei r chemical 
evolu ti on exist ( e - £- IPei & Fall Il995l: iKauffmannl 
I l99d: ljhn"enez et all Il99fit IMo. Mao, fc White! Il999l: 
Pmntzos^^^mssJ^rj_j2000|J 

200ll:lMaller. Prochaska. Somerville. fc Primackl200lLl200l 
Lanfranchi fc FriacafcOOah . Dut they often neglect the effect 
of galaxy mergers, infall of gas, outflows due to supernovae, 
and the interaction with the intergalactic medium (IGM). 
These dynamical processes are expected to have significant 
influence on the chemical evolution of DLAs, and it would 
therefore be ideal to study the chemical evolution of DLAs 
with full cosmological hydro-simulations which treat all of 
these effects dynamically and self-consistently. The earlier 
numerical studies of Gardner et al. (1997a, b; 2001) for 
example did not address the chemical evolution of DLAs in 
detail. 

In our pre vious paper 

(|Nagamine. Springel Sz Hernauistl 120031) . we studied 
the abundance of DLAs using a series of cosmological SPH 
simulations with varying resolution and feedback strength. 
We showed that the ACDM model is able to account for the 
observed abundance of DLAs at redshift z = 3 — 4.5 quite 
well. Another important conclusion of our study was that 
earlier numerical work overestimated the DLA abundance 
significantly due to insufficient numerical resolution, lack of 
efficient feedback processes, and inaccuracies introduced in 
cooling processes when conventional formulations of SPH 
are used, giving rise to an over cooling problem. In our sim- 
ulation methodology, we avoided thes e numerical difficulties 
by u sing a novel version of SPH (|Springel Sz Hernauistl 
2002) that is based on integrating the en tropy as an 
independent th ermodynamic variable (e.g. iLucvl Il977l : 
lHernauistJll99i) . and which takes variations of the SPH 
smoothing lengths self-consistently into account. 

In this paper, we focus on the overall distribution of the 
star formation rate and metallicity of DLAs in the redshift 
range z — — 4.5. We are interested in how these quantities 
appear in our new SPH simulations, in light of our refined 
modelling of star formation and feedback by galactic winds, 
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Run 


Boxsize 


TV 


ra DM 




6 


z end 


wind 


R3 


3.375 


2 x 144 3 


9.29 x 10 5 


1.43 x 10 5 


0.94 


4.00 


strong 


R4 


3.375 


2 x 216 3 


2.75 x 10 5 


4.24 x 10 4 


0.63 


4.00 


strong 


03 


10.00 


2 x 144 3 


2.42 x 10 7 


3.72 x 10 6 


2.78 


2.75 


none 


P3 


10.00 


2 x 144 3 


2.42 x 10 7 


3.72 x 10 6 


2.78 


2.75 


weak 


Q3 


10.00 


2 x 144 3 


2.42 x 10 7 


3.72 x 10 6 


2.78 


2.75 


strong 


Q4 


10.00 


2 x 216 3 


7.16 x 10 6 


1.10 x 10 6 


1.85 


2.75 


strong 


Q5 


10.00 


2 x 324 3 


2.12 x 10 6 


3.26 x 10 5 


1.23 


2.75 


strong 


D4 


33.75 


2 x 216 3 


2.75 x 10 8 


4.24 x 10 7 


6.25 


1.00 


strong 


D5 


33.75 


2 x 324 3 


8.15 x 10 7 


1.26 x 10 7 


4.17 


1.00 


strong 


G4 


100.0 


2 x 216 3 


7.16 x 10 9 


1.10 x 10 9 


12.0 


0.00 


strong 


G5 


100.0 


2 x 324 3 


2.12 x 10 9 


3.26 x 10 8 


8.00 


0.00 


strong 



Table 1. Simulations employed in this study. The box-size is given in units of /i -1 Mpc, N p is the particle number of dark matter and 
gas (hence x 2), raDM and m ga s are the masses of dark matter and gas particles in units of h~ 1 Mo, respectively, e is the comoving 
gravitational softening length in units of /i -1 kpc, and z e nd is the ending redshift of the simulation. The value of e is a measure of spatial 
resolution. From the top to the bottom row, we refer to R3 & R4 collectively as 'R-series', the next 5 runs (03 to Q5) are called 'Q-series', 
D4 fc D5 are called 'D-series ', and G4 & G5 are called 'G-series'. The 'strong- wind' simulations form a subset of the runs analysed by 
ISoringel HernauistJ l2003hft. 



and our new formulation of SPH. In particular, we will ex- 
amine the dependence of our results on feedback strength, 
and also investigate the influence of numerical resolution on 
the physical properties of DLAs. 

The evolution of metallicity in DLAs has previously 
been studied w ith a cosmological hydro-simulation by 
ICen et alJfcOOfl . using an Eulerian code with a fixed mesh- 
size. The methodology of this simulation is very different 
from ours, but it is clearly interesting to compare their re- 
sults to what we find. We discuss this comparison in Sec- 
tion 

The paper is organised as follows. In Section we 
briefly describe the novel features of our simulations, putting 
particular emphasis on the star formation model and the 
treatment of metal production. We show a few examples 
of visual representations of haloes and associated physical 
quantities in Sectional We then discuss the distribution of 
Hi column density as a function of halo mass in Section 0] 
projected stellar mass density in Section |5] projected SFR 
density in Section |H| and projected metallicity of DLAs in 
Section [7| The evolution of the mean metallicity is described 
in Section |S| Finally, we discuss the implications of our work 
m Sectional 



2 SIMULATIONS 

We analyse a large set of cosmological SPH simulations that 
differ in box size, mass resolution and feedback strength, 
as summarised in Table Q] In particular, we consider box 
sizes ranging from 3.375 to 100/i -1 Mpc on a side, with 
particle numbers between 2 x 144 3 and 2 x 324 3 , allow- 
ing us to probe gaseous mass resolutions in the range 
4.2 x 10 4 to 1.1 x 10 9 /i _1 M . These simulations are partly 
taken from a study of the cos mic star formation history by 
ISpringel &; Hernauistl (|2003bl) . supplemented by additional 
runs with weaker or no galactic winds. The joint analysis of 
this series of simulations allows us to significantly broaden 
the range of spatial and mass-scales that we can probe com- 



pared to what is presently attainable within a single simu- 
lation. 

There are three main new features of our simulations. 
First, we use a "conservativ e entropy" formulation of SPH 
(|Springel & Hernauistll20o'3) . This approach has several ad- 
vantages over conventional versions of SPH. Since the en- 
ergy equation is written with the entropy as the indepen- 
dent thermodynamic variable, rather than the thermal en- 
ergy, the 'pdV term is not evaluated explicitly, reducing 
noise from smoothed estimates of e.g. the density. By in- 
cluding terms involving derivatives of the density with re- 
spect to the particle smoothing lengths, this approach ex- 
plicitly conserves entropy (in regions without shocks), as 
well as momentum and energy, even when smoothing lengths 
evolve adaptively, avoiding the problems noted by e.g. Hern- 
quist (1993). Furthermore, this formulation moderates the 
overcool ing, problem present in earlier formulations of SPH 
^(see a lso lYoshida et alJl20()3 : IPearce et al .11 19991 : ICroft et"aD 

Secon d, we employ an effective sub- resolution model de- 
scribed bv lSpringel &: Hernauist (2003a) to treat star forma- 
tion and its regulation by supernovae feedback in the dense 
interstellar medium (ISM). In this model, the highly over- 
dense gas of the ISM is pictured to be a two-phase fluid 
consisting of cold clouds in pressure equilibrium with a hot 
ambient phase. Each gas particle represents a statistical mix- 
ture of these phases. Cold clouds grow by radiative cooling 
out of the hot medium, and this material forms the reser- 
voir of baryons available for star formation. We assume that 
star formation takes place on a characteristic time-scale £*, 
and that a mass fraction f3 of these stars are short-lived and 
instantly die as supernovae so that the local SFR is: 

where p c is the gas density of the cold cloud phase. The 
value of f3 is determined by the initial mass function of the 
stellar population, and here we adopt /3 = 0.1. The star 
formation time-scale £* is taken to be proportional to the 
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local dynamical time of the gas: 

' \ -1/2 



**(p)=*o 



pth y 



(2) 



where the value of to = 2.1 Gyr is chosen to match the Ken- 
nicutt law (1998). Because of its density dependence, the 
star formation time-scale becomes short in very high density 
regions, approximately modelling intense starburst activity. 
In our model, this parameter simultaneously determines a 
threshold density p t h, above which the multiphase structure 
in the gas, and hence star formation, is allowed to develop. 
This physical density is 8.6 x 10 6 /i 2 M©kpc -3 for all simu- 
lations in this study, corresponding t o a comoving baryonic 
overden sity of 7.7 x 10 5 at z = 0. See lSoringel &; Hernauistl 
(2003a) for a description for how p t h is determined self- 
consistently within the model. 

In our simulation code, the mass of a new star particle 
is fixed to ra* = mo/N g , where mo is the initial mass of 
each gas particle, and N g (which we take to be 2 here) is 
the number of "generations" of stars each gas particle may 
spawn. New star particles are created probabilistically with 
expectation value consistent with the current SFR of each 
gas particle. 

Once star formation occurs, the resulting supernova ex- 
plosions are assumed to deposit energy into the hot gas with 
an average of csn = 4 x 10 48 ergs Mq 1 for each solar mass of 
stars formed, or 0.4% of the canonical SN energy of 10 51 ergs. 
The efficiency of SN energy feedback is often described as 
the ratio of the returned SN energy to the rest mass energy 
of the star formed, in which case our SN feedback efficiency 
translates to 6sn/M©c 2 = 2.2 x 10 -6 . The heating rate due 
to SN is 



dt 



(phU h ) 



dp* 
dt 



(3) 



where usn = (1 — 0)f3~ 1 esN is the 'SN energy scale' cor- 
responding to a 'SN temperature' of Tsn = 2/jusN/Sk ~ 
10 8 K, and ph and uh are the density and thermal energy of 
the hot phase of the gas, respectively. 

We assume that SN explosions also evaporate the cold 
clouds, and transfer cold gas back into the ambient hot 
phase. This can be described as 



dpc 
dt 



(4) 



The efficiency parameter A has the density dependence 



-4/5 



(5) 



following iMcKee Ostrikerl (|l977h . The value of A is fixed 
to 1000 by requiring Tsn/^4o = 10 5 K. This evaporation pro- 
cess of cold clouds establishes a tight self-regulation mecha- 
nism for star formation in the ISM, where the ambient hot 
medium quickly evolves towards an equilibrium temperature 
equal to 



u h - 



USN 



+ U c 



(6) 



A + l 

in a time-scale shorter than the star formation time-scale. 
Therefore, it is a good approximation to assume that the 
condition of self-regulated star formation always holds. We 
can then avoid treating the mass exchange between the hot 



and cold phases explicitly. This simplifies the code, makes 
it faster, and reduces its memory consumption. In this 'ef- 
fective' model, which is adopted for all the runs used in the 
present study, the thermal energy of the hot phase is simply 
given by equation jjj, and that of the cold phase is fixed 
to the equivalent of T c = 1000 K. Our result does not de- 
pend on the choice of T c as long as T c <C 10 4 K, because the 
current code does not include molecular hydrogen or metal 
line cooling which would act at temperatures below 10 4 K. 
The mass fraction in the cold clouds for the particles in the 
multiphase medium can be computed by 



Pc 1 . 1 

p 2y 



1 1 

y 4y 2 ' 



where 



y 



t*A ne t(p, Uh) 



p[(3usn - (1 - f3)u c ] ' 



(7) 



(8) 



Here, A ne t(p, u) is the usual cooling function. 

If the physical density of a gas particle is lower than 
pth, the particle represents ordinary gas in a single phase. 
In this case, the neutral hydrogen mass of the particle can 
be computed as 



(p < Pth), 



(9) 



where Xh = 0.76 is the primordial mass fraction of hydro- 
gen, and iVh is the number density of neutral hydrogen atoms 
in units of the total number density of hydrogen nuclei. The 
quantity Nh is computed by solving the ionisation balance 
as a function of density and current UV background flux. If 
the gas density is higher than p t h, the gas particle represents 
multiphase gas, such that the mass of neutral hydrogen can 
be computed as 



mm 



(p > Pth), 



(10) 



where x is the mass fraction of cold clouds given by equa- 
tion (J. 

The simulation also keeps track of metal enrichment, 
and the dynamical transport of metals by the motion of the 
gas particles. Metals are produced by stars and returned 
into the gas by SN explosions. The mass of metals returned 
is AM Z = t/,AM„ where y* = 0.02 is the yield, and AM* 
gives the mass of newly formed long-lived stars. Assuming 
that each gas element behaves as a closed box locally, with 
metals being instantaneously mixed between cold clouds and 
the ambient hot gas, the metallicity Z — Mz/M g of a star- 
forming gas particle increases during one timestep At by 



AZ=(l-%a 



At 



(11) 



When a star particle is generated, it simply inherits the 
metallicity of its parent gas particle. 

Finally, the third important feature of our simulations 
is a phenomenological model for galactic winds, which we 
introduced in order to study the effect of outflows on DLAs, 
galaxies, and the intergalactic medium (IGM). In this model, 
gas particles are stochastically driven out of the dense star- 
forming medium by assigning an extra momentum in ran- 
dom directions, with a rate and magnitude chosen to repro- 
duce mass-loads and wind-speeds similar to those observed. 
The wind mass-loss rate Mw is assumed to be proportional 
to the star formation rate, and the wind carries a fixed frac- 
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DLA 




SFR 





Z 



Figure 1. Projected spatial distribution of various quantities (all in log scale) for the most massive halo of mass Mh a i = 1.7x 10 12 h~ 1 Mq 
at z = 3 in the 'Q5'-run. From top left to bottom right: Nhi, DLAs (log Nhi > 20.3), stellar surface mass density, SFR surface density, 
metal mass surface density, gas metallicity. The size of each panel is comoving zb457/z — 1 kpc from the centre of the halo. The inset in the 
top left panel shows the probability distribution function of lines-of-sight (dn/d log Nhi) for this halo as a function of log Nhi- One can 
see that the majority of the red region is occupied by lines of sight with 11 < log Nhi < 15. 
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tion x of SN energy: 

Mw = r\M* , (12) 
and 

^M w vl = xe SN M*. (13) 

We adopted a fixed value of r\ — 2, and explored two values 
of x 5 a 'weak wind model' with % = 0.25, and a 'strong 
wind model' with x = 1.0 . Solving for the wind velocity 
from the above two equations, the wind models correspond 
to speeds of v w = 242 kins" 1 and 484 kms -1 , respectively. 
This wind energy is not included in csn discussed above, so 
in our simulations (1 + x)csn is the total SN energy returned 
into the gas. 

Most of our simulations employ the 'strong' wind model, 
but for the 10/i _1 Mpc box (runs 03, P3, Q3, Q4, Q5; col- 
lectively called 'Q-series') we also varied the wind strength. 
Therefore, this Q-series can be used to study both the effect 
of numerical resolution and the consequences of feedback 
from galactic winds. The runs in the other simulation se- 
ries then allow the extension of the strong wind results to 
smaller scales ('R-Series'), or to larger box-sizes and hence 
lower redshift ('D-' and 'G-Series'). Our naming convention 
is such that runs of the same model (box-size and included 
physics) are designated with the same letter, with an addi- 
tional number specifying the particle resolution. 

Our calculations include a uniform UV back g round 
radiation field of a modified lHaardt &; Madaul ([l996) 
spec trum, where rei onisation takes place at z ~ 6 
(see |P av^taj2|l999), as suggested by observations (e.g. 
LBecker et alJ 12001') and radiative transfer calculations of 
the impact of the stellar sources in our simulations on 
the IGM (e.g. ISokasian et al.ll2003h . The radia tive cooling 
and h eating rate is computed as described by iKatz et alJ 
(|l996F) . The adopted cosmological parameters of all runs 
are (On, Ov, n b , cr 8 , h) = (0.3,0.7,0.04,0.9,0.7). The sim- 
ulations were performed on the Athlon-MP cluster at 
the Center for Parallel Astrophysical Computing (CPAC) 
at the Harvard-Smithsonian Center for Astrophysics, us- 
ing, a modified version of the parallel GADGET code 
(Springel, Yoshida & WhitejooH) . 



3 VISUAL REPRESENTATIONS 

In Figures [QE1 and El we show the spatial distribution of the 
following quantities for 3 representative haloes in our high- 
est resolution simulation at z = 3 ('Q5'-run): neutral hydro- 
gen column density Nhi, DLAs (simply setting a threshold 
of log Nhi > 20.3 to the total Nhi), projected stellar mass 
surface density, projected star formation rate surface den- 
sity, projected metal mass surface density, and projected 
metallicity. The size and the mass of each halo is given in 
the caption. The pixel size of the picture is taken to be 
the gravitational softening length of the simulation, which 
is 1.23/i _1 kpc for the 'Q5'-run. We defer a discussion of 
quantitative results to the following sections and here sim- 
ply summarise notable features of each figure. 

Figure Q This halo is the most massive one at z — 3 
(Mhaio = 1.7 x 10 12 /i _1 M©), which presumably evolves 
into a group of galaxies at lower redshift. Note that this 
is a very rare halo at z = 3. With our limited boxsize of 



Lbox = 10/i -1 Mpc, only a few such massive haloes (which 
have roughly the space density of Lyman-break galaxies) 
can exist in our simulation box. The neutral hydrogen gas 
distribution in this halo is very extended, and the DLAs are 
widely distributed within the halo, but with a trend of higher 
concentration towards the centre. Obviously, the DLA cross- 
section of this halo is v ery high compared to lower mass 
haloe s (see Figure 2 of iNagamine. Springel. Sz Hernauistl 
2003). The inset in the top left panel shows the proba- 
bility distribution function of sight-lines (dn/d log Nhi) for 
this halo as a function of log Nhi- One can see that the 
majority of the red region is occupied by sight-lines with 
11 < log Nhi < 15, which could be observed in the Ly-a 
forest. As the stellar mass density picture shows, a few high 
concentrations of stellar mass at the centre would probably 
correspond to LBGs at z = 3 (a detailed analyses of LBGs 
in our simulation and their relation to DLAs will be studied 
elsewhere). The projected SFR surface density has a similar 
distribution as the stellar mass, but they do not necessarily 
coincide with each other. The metallicity distribution is still 
more complex. While metals are widely distributed in gen- 
eral, they are also clearly concentrated in the region of high 
SFR. A clear feature is the high metallicity region (white) 
at the centre of the halo, where Nhi is also very high. An 
interesting filamentary structure with low metallicity can be 
seen in the upper left side of the halo, and there is also a 
corresponding feature in the Nhi panel. In this region, there 
is a lot of neutral hydrogen, but the gas has not yet been 
significantly polluted by star formation, as evidenced by the 
low stellar mass density in the same region, explaining the 
low metallicity. 

Figure El This halo of mass M ha i = 2.6 x lO lo /i _1 M 
belongs to a much more abundant class of objects than the 
one shown in Figure Q As can be seen in the panel with the 
stellar mass density (upper middle column), there are 5 or 
6 galaxies in this halo, and the left-most clump seems to be 
going through a merger. The distribution of DLAs, SFR, and 
metal mass all correspond to the five stellar density peaks 
very well. Again, the metallicity is high in these galaxies. 

£ igure 01 This halo of mass M ha i = 2.4 x lO^M© 
is even less massive than the one shown in Figure El by 
an order of magnitude. Two peaks can be seen in the stel- 
lar mass surface density distribution, but only one of them 
(right one) is a DLA galaxy, and only this galaxy is form- 
ing stars. The one on the left presumably converted all of 
its high density cold neutral gas into stars, and therefore 
no DLA and no star formation activity can be seen. The 
metal mass is concentrated in two galaxies, and the metal- 
licity is high in the centre of the halo. The panel on the very 
left shows the probability distribution function of sight-lines 
(dn/dlogNm) through this halo as a function of log Nhi, 
and compares it to the one shown in Figure El One can see 
that the majority of the orange region is occupied by sight- 
lines with 11 < log Nhi < 16, which corresponds to the 
column densities of the Ly-a forest. 



4 Hi COLUMN DENSITY VS. HALO MASS 

In Figure 31 we show the Hi column density of DLA 
sight-lines against the dark matter halo mass in which 
DLAs reside. The calculation of the Hi column den- 
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Figure 2. Projected spatial distribution of various quantities for a halo of mass M^ab = 2.6 x lO lo /z _1 M0 at z = 3 in the 'Q5'-run. 
Left column: Nhi (top) and DLAs (bottom). Middle column: stellar surface mass density (top) and SFR surface density (bottom). Right 
column: metal mass surface density (top) and gas metallicity (bottom). The size of each panel is comoving ±112/i — 1 kpc from the centre 
of the halo. 
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Figure 3. Projected spatial distribution of various quantities for a halo of mass Mhaio = 2.4 x 10 9 h~ 1 MQ at z = 3 in the 'Q5'-run. 
The size of each panel is comoving zb57/i _1 kpc from the centre of the halo. The graph on the top left shows the probability distribution 
function of lines-of-sight (dn/d log Nhi) for this halo and the one shown in Figure 1^1 as a function of log Nhi- One can see that the 
majority of the orange region is occupied by lines-of-sight with 11 < log Nhi < 16, which could be observed in the Ly-cn forest. The other 
panels are ordered in the same way as in Figure El 
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Figure 4. Hi column density vs. dark matter halo mass at z = 3. Each point in the figure represents one line-of-sight. Contours are 
equally spaced on a logarithmic scale and the lowest contour contains 95% of all the points. The highest resolution run in this figure is 
Q5 (bottom left panel). 



sity is carried out exactly the s ame w ay as we did in 
Nae amine. Spring el. Hernanistl ((2003), which we sum- 
marise briefly here. First, we identify dark matter haloes by 
applying a conventional friends-of-friends algorithm to the 
dark matter particles in each simulation. After dark matter 
haloes are identified, we associate each gas and star parti- 
cle with their nearest dark matter particle, including them 
in the particle list of the corresponding haloes, when ap- 
propriate. Then, for each halo, a uniform grid covering the 
entire halo, with grid-size equal to the gravitational soften- 
ing length, is placed at the centre-of-mass of the halo. The 
neutral mass of each gas particle is smoothed over a spheri- 
cal region of grid-cells, weighted by the SPH kernel. We then 
project the neutral gas in the halo onto a plane, and obtain 
the column density of each grid-cell in this plane. 

We regard each grid-cell on the x?/-plane as one sight- 
line, as we have no power to resolve smaller scales than the 
gravitational softening length of each simulation. Each point 
in the figure represents one sight-line. Contours are in equal 
logarithmic scale, and the lowest contour contains 95% of 
all the points. 

A feature to note here is that a wide range of Nhi values 
correspond to massive haloes at z — 3. This means that the 



massive haloes that host LBGs at z — 3 also host a wide va- 
riety of DLAs. Note that, since we identified the dark matter 
haloes with a simple friends-of-friends algorithm, subhaloes 
in massive haloes are not separated from the parent halo 
in this figure. There is also a weak trend that less massive 
haloes in general host lower Nhi systems, which is naturally 
expected. Low mass haloes with Mhaio < 1O 9 M tend to 
host only one DLA with lower Nhi at the centre of the halo. 



5 PROJECTED STELLAR MASS DENSITY 
AND DLAS 

5.1 Redshift z = 3 

In Figure |5] we show the distribution of projected stellar 
mass density (in proper units) of DLAs as a function of 
Hi column density. Each point in the figure represents one 
line-of-sight. Contours are in equal logarithmic scale, and 
the lowest contour contains 95% of all the points. Sight-lines 
with no stars are shown as crosses at the bottom of each 
panel. As a reference, the solid horizontal line indicates an 
estimate of surface mass density of stars in the Milky Way, 
Sstar = 50 M pc -2 , although this number is highly uncer- 
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Figure 5. Projected stellar mass density vs. Hi column density at z = 3 (both in proper units, not comoving; same for the rest of the 
figures). Each point in the figure represents one line-of-sight. Contours are equally spaced on a logarithmic scale and the lowest contour 
contains 95% of all the points. The highest resolution run in this figure is Q5 (bottom left panel). The long-dashed line indicates the 
relation of the projected star-to-gas mass ratio ^star/^pNhi = 3. See text for a discussion of the other 2 horizontal lines in the figure. 
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Figure 6. Projected stellar mass density vs. Hi column density at z = 4.5. The long-dashed line indicates the projected star-to-gas mass 
ratio Sstar/^pNhi = ^ (instead of 3 in Figure 0. The other two horizontal lines are the same as in Figure |3 
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Figure 7. Projected stellar mass density vs. Hi column density 
at z = 1 and z = 0. The long-dashed line indicates the projected 
star-to-gas mass ratio Sstar/^pNhi = 10 at z = 1 (top 4 pan- 
els), and 20 at z = (bottom 2 panels). The other two horizontal 
lines are the same as in Figure 151 



tain because of the difficulty in constraining the contribution 
of da r k matter in the disc (iB ahcall 1984; Kuiiken^^^ilmore 
Il98l Il99ll: iBahcall. Flvnn. fc Gould Il992l : iFlvnn fc Fuchs 
119941 : iBinnev fc Merrifiekl ll998l Most recent studies using 
the Hipparcos satellite data suggest that there is no com- 
pelling e vidence for significant amounts of dark mat ter in 
the disc dCreze et al.lll998l : iHohnberg Flvnnl l200(tl . The 
short- dashed horizontal line indicates a reference value for 
5^ star if stars of total mass 10 n M Q are distributed in a cir- 
cle of radius 10 kpc. As a reference, we show the projected 
star-to-gas mass ratio of Sstar/ ^pNni = 3 with the long- 
dashed line, where m p is the proton mass. The points are 
distributed on both sides of this line at 20 < logNm < 22.5, 
but they preferentially lie above this line at logNm > 22.5 
in all of the panels, suggesting efficient star formation at 
these high column densities. Stripes seen at low values of 
Sstar show the discreteness due to the limited stellar mass 
resolution of stars in the simulation. 

Comparing the top 3 panels of Figure|5| '03' (no wind), 
'P3' (weak wind), and 'Q3' (strong wind), gives us an idea 
about the effect of the feedback by galactic winds. The runs 
with no winds (03) and weak winds (P3) contain more sight- 
lines with high Sstar when compared to 'Q3' (strong wind 
run, upper right corner of the panel). This is because the gas 



in '03' can cool and form stars continuously without being 
blown out by galactic winds. 

As we increase the numerical resolution from Q3, Q4, 
and to Q5 (all strong wind runs), the number of sight-lines 
with log Nhi > 22 increases, because higher resolution runs 
can resolve higher gas density. Otherwise, the distribution of 
points in Q3, Q4, and Q5 looks similar. Comparison between 
D4 & D5, and G4 & G5 also shows a similar trend. 

As the boxsizes become larger from the Q-Series (runs 
with 10/i _1 Mpc box) to the D- and G-series, the bottom 
end of the contour rises due to degrading mass resolution. A 
larger simulation box (G5) contains more high-mass haloes 
than a smaller box (Q5); therefore G5 contains more sight- 
lines with high Sstar values. 

5.2 Redshift z = 4.5 

In Figure El we show the same quantities as in Figure |5] but 
at redshift z = 4.5. Horizontal solid and short-dashed lines 
are the same as in Figure The long- dashed line here shows 
the projected star-to-gas ratio of Sstar/^pNhi = 1. The 
star-to-gas mass ratio at z — 4.5 is lower than that at z — 3 
by about a factor of a few. 

The highest resolution run in this figure is R4 (bot- 
tom right panel), but the boxsize of this run is small 
(3.375/i -1 Mpc). A more reliable distribution can be ob- 
served in the Q5 run. 

5.3 Lower redshift 

In Figure Q we show the same quantities as in Figure |5| but 
at redshift s z — 1 and z — 0. The two horizontal lines are 
the same as those in Figure|5] The long-dashed line indicates 
the projected star-to-gas mass ratio SsTAR/mpNni = 10 at 
z — 1 (top 4 panels), and 20 at z — (bottom 2 panels). 

The observational estimate of the local star-to-gas 
mass ratio in the Milky Way is Estar/wIpNhi ~ 10. 
This is obtained by dividing the estimate of Sstar ~ 
50 Mqpc -2 (assuming that there is no s ignifican t con- 
tribution of dark matte r in t he disc) (pahcall| [l984|: 
Kuiiken & Gilmore [l989l. Il99ll: IBahcall. Flvn^ fc Gould 



19921: IFlvnn fc Fuchsl Il994t IBinnev fc Merrifieidl IK 

Creze et al.lll998l : lHolmberg fc Flvnnl200(t ) by the co ntribu- 
tion f rom the interstellar medium, Eism ~ 5 M^pc -2 (jPamd 
1993). Given the uncertainties in these numbers, the values 
we find in the simulation are consistent with the observed 
range. 

An interesting evolutionary feature in the distribution 
can be observed as a function of redshift. At z = 3, the 
distribution is wide in the horizontal direction, but as the 
redshift decreases, the distribution becomes narrower and 
elongated in the vertical direction, and the points move to 
higher Sstar and lower Nhi values as the neutral gas is 
converted into stars inside the galaxies. 



6 STAR FORMATION RATE OF DLA 
GALAXIES 

In Figure |H] we show the projected star formation rate sur- 
face density Ssfr (in proper units of M yr -1 kpc -2 ) as a 
function of Hi column density. The SFR of each gas particle 
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Figure 8. Projected star formation rate density vs. Hi column density at z = 0, 1, 3, Sz 4.5. Each point in the figure represents one 
line-of-sight. Contours are equally spaced on a logarithmic scale and the lowest contour contains 95% of all the points. The shaded region 
in the Q5 panel roughly indicates the region of recent observational data by Wolfe et al. (2003b). 



is smoothed over a spherical region of grid-cells, weighted 
by the SPH kernel. This ensures that rate of star formation 
is conserved when added up along the ^-direction to obtain 
the projected SFR density. 

Each point in the figure represents one line-of-sight. 
Contours are equally spaced on a logarithmic scale, and the 
lowest contour contains 95% of all the points. The distri- 
bution was truncated at logNm = 20, so the cutoff of the 
contours at logNm ~ 20 is artificial. The three solid lines 
in each panel show the allowed region by the following em- 
pirical relation between the projected SFR densit y and the 
Hi co lumn density, known as the Kennicutt law (|Kennicuttl 

W 



S SFR = (2.5 ±0.7) x 10" 4 (- 



Nb 



25 x 10 20 cm- 2 , 

M©yr _1 kpc~ 



(14) 
(15) 



The projected SFR of DLAs in the simulation agrees 
well with the Kennicutt law. It may at first be surprising 
that the simulated result at z = 3 follows the Kennicutt 
law so well, which empirically is only established for z = 0. 
We will discuss this point further in Section There is a 
slight trend in all panels for the distribution to deviate to 
higher values of Ssfr at 20 < logNm < 21 relative to the 
Kennicutt law, but the reason for this deviation is not clear. 

The shaded area in the panel of Q5 in Figure |H| 
roughly indicates the region of recent observational data by 
Wolfe, Prochaska, & Gawiser (2003b). Their data points are 
close to the simulated result and roughly falls on the Kenni- 
cutt law, but the individual measurements do not with some 



scatter. Wolfe et al. (2003b) argue that this scatter is not 
surprising because their measurement of projected SFR is 
an average over a scale of a couple of kpc over which Nhi is 
expected to fluctuate. 

As the wind strength increases from 03 to Q3, the 
number of high Nhi points decreases for the same reason 
described in Section 15.11 but otherwise the distribution re- 
mains very similar for different wind models. This is because 
even though the Hi content of galaxies can be changed by a 
galactic wind, the star formation recipe in the code gives the 
same SFR for a given Hi content. Varying the wind strength 
therefore only moves the points along the Kennicutt law. 

Increasing the numerical resolution from Q3, Q4, and 
then to Q5, slightly increases the number of sight-lines at 
high Nhi values, but the distribution is basically identical. 
In the G-series, one sees a sprinkle of points that fall away 
from the Kennicutt law at around 20 < log Nhi < 21, which 
is perhaps due to a lack of resolution in these large boxsize 
runs. 

At z = 4.5 (only Q5 is shown), the distribution of points 
seems to be tighter than that at z = 3. An interesting fea- 
ture seen at z — 4.5 is the significant number of points falling 
away from the Kennicutt law at log Nhi ~ 20.5. This is in- 
dicating that there are many more DLA sight-lines with in- 
efficient star formation at z — 4.5 compared to z = 3, which 
is presumably related to the sudden decline in the DLA 
cross-section at around Mhaio = 10 8 ' 3 h ~ 1 M<7 ) , as found in 
iNagamine. Springe!. Sz Hernauistl (|2p03). At z = 4.5, many 
low-mass haloes still have not accumulated enough gas for 
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Figure 9. Gas metallicity vs. Hi column density at z = 3. Each point in the figure represents one line-of-sight. Contours are equally 
spaced on a logarithmic scale. Current observational data points fall into the shaded region that is drawn in the bottom left panel (Q5). 



them to radiatively cool, and the neutral gas content and 
the SFR is significantly lower than in those at z = 3. 

At low redshift (z = & 1), the distribution starts to 
shrink towards lower Ssfr and lower Nhi values. This can 
be largely understood in terms of the global evolution of the 
star formation efficiency in haloes, and its variation with 
haloes mass. As the redshift decreases, a larger fraction of 
the total star formation rat e in the universe is contribute d by 
ever more massive haloes (Spri ngel fc Hernqu ist 2003b|). In 
addition, at lower redshifts, haloes are less dense, and their 
gas has been heated by prior star form ation, such that cool - 
ing becomes comparatively inefficient ([Blanton et al.l ll999). 
leading to a reduction in the star formation rates. 



7 GAS METALLICITY OF DLAS 
7.1 Redshift z=3 

In Figure E| we show the projected gas metallicity as a func- 
tion of Hi column density at z = 3. The projected metallicity 
was calculated as follows: First, the metal mass and gas mass 
of each gas particle was smoothed over a spherical region of 
grid-cells, weighted by the SPH kernel. We then added up 
gas and metal mass independently along the z-direction, and 
finally divided the sum of the metal mass by the sum of the 
gas mass for each sight-line. We have also tried the same 
calculation by weighting with the neutral hydrogen mass in- 
stead of using the total gas mass, but the result was almost 
identical. We also confirmed that the result is robust even 
when we select out only the gas particles in high density 



clumps, therefore the results are not affected by the diffuse 
gas in the halo when projected along the ^-direction. 

Each point in the figure represents one sight-line, and 
the contours are equally spaced on a logarithmic scale. The 
solid square symbols give the median value in each log Nhi 
bin, with error bars indicating the quartiles on both sides. In 
all panels, the median metallicity increases as Nhi increases, 
reaching solar metallicity at log Nhi ~ 23. This trend is 
expected because star formation is more vigorous in high 
Nhi systems, as we saw in Section |S] 

In the '03' run (no wind), the contours extend to higher 
metallicity than in the 'Q3' run (strong wind). This is be- 
cause in 03 star formation is not suppressed by a galac- 
tic wind. Therefore, the gas is more enriched with metals 
due to a higher average star formation rate. The values 
of the median metallicity at intermediate column densities 
(20 < log Nhi < 22) are similar for all runs in the Q-series 
(03, P3, Q3, Q4, and Q5). 

Increasing the resolution from Q3, to Q4, and finally 
to Q5, increases the number of points with high metallicity 
at high Nhi values, but otherwise the overall distribution is 
very similar. 

It is clear that simulations with larger boxsizes, such as 
those of the D- and G-series, underestimate the metallicity 
for systems with log Nhi < 22 at z = 3 due to lack of res- 
olution. This is consistent with the finding that in the runs 
of the D- and G-series the Hi column density distribution 
function at z — 3 is not reliably estimat ed, as was shown in 
iNagamine. Springel. &; H ernquist! (|2QQ3h . 

An important result is that the median metallicity we 
find for DLAs in our simulations is much higher than the ob- 
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Figure 10. Projected SFR vs. gas metallicity at z = 3 for the 
'Q5'-run. Each point in the figure represents one line-of-sight. 
Contours are equally spaced on a logarithmic scale. The shaded 
area indicates the region of observed data points by Wolfe et 
al. (2003b). 



served value. Note that observers typically find DLAs with 
20 < log Nht < 22 t o have a meta l licity o f Z/Z^ ~ 1/30 (e.g. 
[ Boisse et all llQQfil IPettini et all llQQflt IProchaska Wolfe! 
2000), as indicated by the shaded region in the panel of 
'Q5' (bottom left). We will discuss the implications of this 
discrepancy further in Section 

In Figure HUl we show the distribution of SFR as a func- 
tion of metallicity at z = 3 for the 'Q5'-run. The recent 
observa tional data points by IWolfe. Prochaska. Gawiserl 
(2003b) fall into the shaded area. The shape of the simu- 
lated distribution is easy to understand: Because Esfr is 
tightly correlated with Nhi, as we saw in Figure |H] (the Ken- 
nicutt law) , the distribution seen in Figure 1101 is simply a 
reflection of Figure around a diagonal line. 



7.2 Redshift 



4.5 



In Figure 1111 we show the same quantities as in Figure |5| 
but for z — 4.5. A noteworthy feature is that the median 
metallicity at z = 4.5 is lower than that at z = 3 by a 
factor of 2 — 3 at column densities log Nhi < 22. This is 
consistent with the fact that the star-to-gas mass ratio at 
z = 4.5 is lower than that at z = 3 by a similar factor. 
The trends observed for increasing numerical resolution and 
wind strength are similar to those at z — 3. 



7.3 Lower redshift 

In Figure we again analyse the same quantities as in 
Figure but at redshift s z = 1 and z = 0. The most reliable 
of our runs at z = 1 is 'D5', which shows a median metallicity 
that is higher by a factor of ^ 2 compared to that at z = 3. 
A factor of 3 increase is observed for 'G5'. 



We note that the evolution in the median metallicity 
seen from z — 1 to z — is quite uncertain. In fact, for 
the 'G5' run, the 3 median points at 20 < log Nhi < 21 go 
down by a factor of 2 from z = 1 to z = 0. This is because 
in a simulation with a large boxsize, such as 'G5', the mass 
resolution is low and many haloes are just forming at low 
redshift. Therefore, the number of points at 20 < log Nhi < 
21 with low metallicity is still increasing from z = 1 to z = 
0, as evidenced by the broadening of the contours towards 
lower metallicity. This is consistent with the fact that the 
value of Qhi in 'G5' slight ly increases from £ = ltoz = 0, as 
was sh own in Figure 1 of iNagamine. Springel. &; Hernauistl 
(2003). 



8 MEAN METALLICITY EVOLUTION 

In FigurelTSI we compare observations of the evolution of the 
Nm-weighted mean metallicity of DLAs with results mea- 
sured from our simulations. Here, by "Nm-weighted" , we 
mean summing up all the metal mass in DLA sight-lines 
and dividing it by the total gas mass in all DLA sight-lines. 
The symbols connected by the solid lines are the direct re- 
sults of our simulations for DLAs in the Q5, D5, and G5- 
runs. We have also tried to do the calculation in a different 
manner; simply taking the mean of the metallicities of all 
lines-of-sight in the simulations; i.e. (Z/Zq) = j E(Z/Zq)/N, 
where N is the total number of lines-of-sight. This method 
gives a mean metallicity which is lower by ~ 0.1 dex than 
the one obtained by the former method. This is perhaps be- 
cause of the wide range of values that DLA metallicities can 
take, as seen in Figure |5| Our conclusion regarding the DLA 
metallicity does not change depending on which calculation 
method is used. 

Column density-weig ht ed m ean observational data 
points bv | Pet tini_et all (|l999l) (open squares) and 
IProchaska et alJ (I2003T) (crosses ) are a lso shown. Note that 
the err or bars of IPet tini et all (|l999F) points are 1-cr, and 
those of lProchaska et al.l (|2003i) are 95% confidence level un- 
certainty, both determined from a bootstrap error analysis. 
The long-dashed line is the best-fit line obtaine d by perform- 
ing a least- square fit to the mean points of lProchaska et all 
(2003), which exhibits a mild evolution with a slope of —0.30 
as a function of redsh i ft [log (Z/Z ) -0.30z - 0.44]. Note 
that IProchaska et alJ |2003) performed a least-square fit to 
the central values of the error bars, so they report a slope of 
-0.26 which is slightly different from our value, but the two 
agrees with each other within 1-cr. 

The mean metallicity of DLAs in our simulations is 
again higher than current observations by about a factor 
of 3 — 5. On the other hand, the mean metallicity of the 
entire simulation box (i.e. total metal mass divided by the 
total gas mass) is lower than the observed DLA metallic- 
ity. This makes sense because the box-mean includes low 
metallicity IGM that is not polluted as strongly as galaxies 
themselves. This simple consideration also illustrates that 
the metallicity values of DLAs can be quite sensitive to de- 
tails of how metals are transported and mixed with gas in 
the haloes of galaxies and with gas in the IGM. In fact, we 
find that our simulations predict the correct metallicity of 
Lyman-a forest when the winds are included (Springel et al. 
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Figure 11. Gas metallicity vs. Hi column density at z = 4.5. Each point in the figure represents one line-of-sight. Contours are equally 
spaced on a logarithmic scale. 



2003, in preparation), and that the box-mean metallicity is 
biased higher relative to that of the low-density IGM. 

Residual differences between the results of Q5, D5, and 
G5 are due to numerical resolution. Higher resolution runs 
(like Q5) can resolve smaller haloes at high redshift than 
lower resolution runs (D5, G5). Therefore, Q5 has more star 
formation, and consequently more metal production, at high 
redshifts than D5 and G5. 

Another noteworthy result is that the rate of the 
evolution of mean metallicity in our simulations appears 
to be consistent with the most recent observe d rate by 
IProchaska et alJ (|2QQ3h . In the simulation of ICen et all 
(2002|), a similar agreement was obtained. We will discuss 
this more in the following section. 



9 DISCUSSION 

We have used a series of state-of-the-art hydrodynamic sim- 
ulations to study the distribution of star formation rate and 
metallicity of DLAs. Our simulations suggest that the me- 
dian projected star-to-gas mass ratio of DLAs is about a few 
at z = 3, and then evolves to higher values of c± 10 at z = 1 
and ~ 20 at z = 0. However, we note that our low-redshift 
results should be regarded with caution because they are 
based on simulations with coarser resolution than those at 
high- redshift. Nevertheless it is encouraging that the star- 
to-gas mass ratio we find at z = is close to observational 
estimates for the Milky Way (see Section . 

The projected SFR density as a function of neutral hy- 
drogen column density follows the Kennicutt law well at 



all redshifts. At first sight, this may seem non-trivial and 
perhaps somewhat surprising. However, the star formation 
model adopted in the simulations depends only on local 
physical quantities, e.g. the local gas density. The free pa- 
rameter of the model was chosen to reproduce the Kenni- 
cutt law in isolated disk galaxies at low redshift, but was 
kept fixed as a function of time. We hence implicitly as- 
sumed that the Kennicutt law holds at all redshifts, and our 
simulation results reflect this assumption. 

A recent study bv iKravtsovl ([2003F) also suggests that 
the projected SFR follows the global 'Schmidt' law, and that 
it does not depend on redshift. He argues that the global 
Schmidt-law originates from a generic power-law distribu- 
tion of the high-density tail of the probability distribution 
function of gas density. If the star formation at high-redshift 
follows similar physical process as in the Local Universe, it 
is not unreasonable that the projected SFR should follow 
the Kennicutt law well even at high redshift. 

The Hi column density distribution in our simulations 
extends up to values as high as log Nhi ~ 23, which is in the 
range of Hi column densities seen in the nuclei of starburst 
galaxies in the Local Universe. This suggests that the intense 
star- formation seen in some of the brightest simulated galax- 
ies at z — 3 is similar to the star formation in local starburst 
galaxies. These simulated galaxies should then correspond 
to the observed Ly man-break galaxies, which are often esti- 
mated to have SFRs of ~ 50 — 100M©yr _1 , or even higher. 
The spatial correlation between Lyman-break galaxies and 
DLAs is hence of significant interest to elucidate this con- 
nection further, and we will study this issue in future work. 

While the star-to-gas mass ratio and the projected SFR 
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Figure 12. Gas metallicity vs. Hi column density at z = and 
1. Each point in the figure represents one line-of-sight. Contours 
are equally spaced on a logarithmic scale. 



density at z — 3 in the simulations seems to be plausible, 
the median metallicity of DLAs at z = 3 appears to be too 
high (Z/Zq ~ 1/3) compare d to typically observ ed values of 
Z/Zq ~ 1/30 for DLAs (e.g. IPettini et alJll99<A . There are 
a number of possible explanations for this problem, and we 
will briefly discuss some of the most prominent possibilities. 

Since the SFRs and metallicities of galaxies in the sim- 
ulation at z = 3 seem plausible, one possible reason for 
high metallicity in DLAs is that the feedback by galactic 
winds is not efficient enough in blowing out metals from 
DLAs. Clearly, if the feedback by winds were stronger, then 
star formation and hence metal creation in DLAs would be 
more strongly suppressed. In fact, the SFR of systems with 
20 < logNm < 21 deviates slightly from the Kennicutt law 
to higher values in our simulation. Reducing the SFR of 
these systems would bring them down to values consistent 
with the Kennicutt law. 

However, simply making the winds stronger and blow- 
ing out more gas will not necessarily decrease the metallic- 
ity of DLAs much, because in our current simulation model, 
the wind transports away metals and gas at the same time, 
i.e. the wind's initial metallicity is assumed to be equal to 
that of the gas of the DLA, leaving the ratio of metal and 
gas mass in the DLA unchanged. 

It is however quite plausible that the wind is metal- 
loaded compared to the gas in the DLA, as is for example 
suggested by simulations of SN explosions (e.g. MacLow & 
Ferrara 1999, Bromm et al. 2003). After all, the ejecta of SN 
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Figure 13. Evolution of the mean gas metallicity of DLAs as 
a function of redshift. The symbols connected by solid lines give 
the direct simulation results for the Q5, D5, and G5-runs. The 
short-dashed lines at the bottom of the figure are the mean metal- 
licities of the entire simulati on boxes. Data p oints with error 
bar s show observations by IPettini et alJ fll999T) (open squares) 
and lProchaska et alJ |2003^ (crosses). The long-dashed line is the 
best-fit line obtained by performing a least-square fit to the mean 
points of lProchaska et alJ |2003ft . which exhibits a mild evolution 
with a slope of —0.30 as a function of redshift. 



are heavily enriched and inject large parts of the energy that 
is assumed to ultimately drive the outflow. If the mixing with 
other DLA-gas is not extremely efficient before the outflow 
occurs, it can then be expected that the wind material has 
potentially much higher metallicity than the DLA, thereby 
selectively removing metals. 

A related possibility concerns the metallicities of the 
cold and diffuse phases of the DLA. In the present study, 
we assumed that metals are always efficiently and rapidly 
mixed between the gas of the cold clouds and the ambient 
medium, such that there is a homogeneous metal distribu- 
tion in the DLA (operationally, we used only a single metal- 
licity variable for each gas particle, reflecting this assump- 
tion). However, this assumption may not be fully correct. If 
the metals were preferentially kept in the hot phase of the 
ISM after they are released by SN, then they would not be 
observed in the cold gas that is responsible for the DLAs. 
Since we did not track the metal distribution in cold and 
hot phases separately in the current simulations, we may 
then have overestimated the amount of metals in DLAs by 
counting those in the hot phase as well as those in the cold 
phase. Note that a more detailed tracking of metals in the 
simulation, separately for hot and cold phases of the gas, 
could in principle be done easily on a technical level. The 
difficulty however lies in obtaining a reasonable description 
of the physics that governs the exchange of metals between 
the different phases of the ISM, something that is presently 
not attainable from either observation or theory. 

The viability of the feedback model in the simulations 
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can also be tested by comparing with observations of the 
Lyman-a forest, which is generated by systems of much 
lower column density than the DLAs studied in this work. 
Curiously, an analysis using the current simulation series 
(Springel et a l. 2003 , in preparation), as well as a study by 
iTheuns et alJ (|2Q02T) . suggest that the spectral features of 
the Ly-a forest are not significantly affected by the feed- 
back from galactic outflows, despite the fact that the wind 
strength is taken to be on the 'strong side' in these stud- 
ies, and despite the fact that a non- negligible fraction of the 
IGM volume is heated by the winds. Note that in our simu- 
lations with a strong wind model, the Hi mass density in the 
entire simulation box is s omewhat lower than suggested by 
obser vational estimates (fNagamine. Springel. &; Hernauistl 
2003), therefore it appears problematic to increase the wind 
strength even beyond the present value. 

The high metallicity of DLAs may also be related to 
the steep luminosity function of galaxies in our SPH simu- 
lations. A preliminary analysis using a population synthesis 
model shows that the luminosity function still has a very 
steep slope at the faint end even at low redshift, similar 
to that of the dark matter halo mass function. (It is not 
obvious whether the steep faint end in the simulation at 
z > 1 is a problem, because it is not well constrained ob- 
servationally at z > 1 yet.) This means that the formation 
of low-mass galaxies in our simulation was not suppressed 
enough, or equivalent ly, that star formation was too efficient 
in low-mass haloes. Stronger winds may help to alleviate 
this problem, but it appears unlikely that our present feed- 
back model can solve it satisfactorily simply by adopting a 
higher efficiency parameter for feedback. It is more plausible 
that additional physical processes need to be considered in 
a more faithful way. One simple possibility for this is related 
to the UV background field, which is turned on by hand at 
z — 6 in our present simulations to mimic reionisation of the 
Universe at a time when the first Gunn - Peterson troughs i n 
spectra to distant quasars are observed (Beck eret al.ll2001r) . 
However, it is possible (in fact suggested by the WMAP 
satellite) that the Universe was reionised at much higher 
redshift. The associated photoheating may have then much 
more efficiently impaired the formation of l ow- mass galax- 
ies th an in our present simulations (but see Diik stra et alJ 
2003). We plan to explore this possibility in future work by 
adopting different treatments of the UV background radia- 
tion field. 

It is also interesting to compare our result s with other 
hydrodynamic simulations of DLAs. ICen et all (|2002f) have 
recently studied the metallicity of DLAs in an Eulerian hy- 
drodynamic simulation of boxsize Lbox = 25/i _1 Mpc. The 
median metallicity as a function of Hi column density in 
their simulation is lower than ours by a factor of ~ 5, in 
better agreement with observations (but still higher than 
observations by about a factor of a few) . This difference be- 
tween the two simulations could be due to the difference 
in the efficiency of the SN feedback, which in turn origi- 
nates in the widely different spatial resolutions of the two 
simulations. The Eulerian simulation of Cen et al. used a 
fixed grid, where the energy and the metals released by SN 
were dumped into the gas of those local cells that showed 
star formation activity. Since the comoving cell size of their 
simulation was 32.6/i -1 kpc, metals were thus distributed ef- 
ficiently over this spatial scale by construction, diluting the 



metal mass relative to the gas mass in the cell for a given 
amount of stars formed. On the other hand, in our present 
SPH simulations, the gravitational softening length was typ- 
ically only a few kpc comoving, such that the gas particles 
would in fact have to travel over a large number of resolu- 
tion elements in order to distribute the injected metals over 
similar spatial scales as in the mesh simulation. 

Despite the differences in the metallicity values of DLAs 
measured in the two simulations, there are also important 
common features. One of them is the existence of high Nhi 
systems with high metallicity, which are absent in current 
observations. The existence of such a population of DLAs 
may pe rhaps be a gener ic prediction of current CDM simu- 
lations. |Ce^^^l] (|2002fl invoked a dust obscuration effect in 
order to reconcile their result with observations, and argued 
that their column density distribution, Hi mass density, and 
D LA- metallicity would all become consistent with observa- 
tions, provided a sizable dust extinction is assumed. How- 
ever, it is pre sently not very clear how strong dust effect s 
really are (e.g.Sj et al.l200ll:IProchaska WolfeEoQ^. 
and the solution could rather lie in a more adequ ate treat- 
ment o f star formation and super nova feedback, as lCen et alJ 
( 2002) also warn. For example, ISchavd (2001) argues that 
the conversion of neutral hydrogen atoms into a molecular 
form, which we have not yet implemented in our simulations, 
would introduce a physical limit to the highest Nhi that 
DLAs can attain. This process would eliminate the highest 
Nhi systems, but would not reduce the metallicity of low 
Nhi systems because it would make the star formation even 
more efficient than it is in current simulations. If future cos- 
mological simulations with a more sophisticated modelling 
of star formation and SN feedback confirm the existence of 
high- metallicity, high-Nm systems, then they could turn into 
an interesting challenge for the CDM model. 

Another important common feature found in both sim- 
ulation studies is the mild evolution of the mean metallicity 
of DLAs with redshift. It is encouraging that the two entirely 
different simulation methodologies agree on the rate of this 
evolution, which in turn is consistent with the observed rate. 
Note that processes such as galaxy mergers, gas infall, and 
outflows all play an important role in the chemical evolution 
of DLAs, and these processes, which often cannot be treated 
accurately in semi-analytic models of the chemical evolution 
of DLAs, are included dynamically and self-consistently in 
both simulations. 

In c onclusion, together with our previous re sults pre- 
sented in lNagamine. Springel. Sz Hernauistl (|2003), we have 
shown that the DLAs found in our simulation series have 
many plausible properties. In particular, they are in good 
agreement with recent observations of the total neutral 
Hi mass density, the Hi column density distribution func- 
tion, the abundance of DLAs, and the distribution of SFR 
in DLAs. However, our simulated DLAs show typically con- 
siderably higher metallicity than what is presently observed 
for the bulk of these systems. This likely indicates that 
metal transport and mixing processes have not been effi- 
cient enough in our simulations. It will be interesting to 
study more sophisticated metal enrichment models in future 
simulations in order to further improve our understanding 
of the nature of DLAs in hierarchical CDM models. 
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